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SUMMARY 


A previous analysis of the creep behavior of a slightly curved pin- 
ended H-section column under constant load is extended to the slightly 
curved solid rectangular- section column. The analysis leads to a dif- 
ferential equation for the plastic strains at the midheight cross 
section. The form of the equation indicates the significant parameters 
which may be useful in plotting test data on the creep life of columns. 
These are a lifetime parameter t' cr , an initial- straightness param- 
eter S or S', and the ratio of the average applied stress to the Euler 
stress a/cr-g. A numerical method of solving the differential equation, 
suitable for use with a high-speed digital computer, is described, and 
typical computed results are given. The existence of a finite lifetime, 
although not evident from the differential equation, is argued intuitively 
and confirmed by the numerical computations . 


INTRODUCTION 


The high temperatures that can develop in aircraft during supersonic 
flight make it important to consider the possible limitations due to 
creep on the useful service life of the structural components. In a 
previous paper (ref. 1), the creep of a slightly curved pin-ended ideal- 
ized H-section column under a constant load and constant temperature was 
studied theoretically. The material of the column was characterized by 
an assumed creep law, for constant uniaxial canpressive stress and con- 
stant temperature, of the form 


e 


E 


+ Ae 


Bcr t K 


where e is the total compressive strain, a is the constant compressive 
stress, t is the time after application of the stress, and E, A, B, 
and K are material constants. This form was selected because it 
applied to at least two alloys: the 75S-T6 aluminum alloy at 600° F 

(ref. 2) and a low-alloy steel at 800° F and, possibly, 1,100° F (ref. 3). 
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Shanley's engineering hypotheses of creep (ref. L) were used in order to 
generalize the assumed constant-stress creep law to cover situations in 
which the stress varies with time, the condition encountered in the 
fibers of a column undergoing continuous lateral deflection because of 
creep . 

In the present paper, the analysis of reference 1 is extended to 
the solid rectangular- section column. Except for the shape of the cross 
section, the assumptions of reference 1 are retained. In addition to 
those cited previously concerning material behavior, these assumptions 
are: (a) that the initial shape is half a sine wave, (b) that the load 

on the column is applied rapidly enough so that negligible creep occurs 
during the loading period (that is, elastic behavior on loading) but 
not so rapidly that dynamic effects have to be considered, (c) that the 
shape of the deforming column remains at all times sinusoidal, (d) that 
plane sections perpendicular to the (curved) center line of the column 
before loading remain plane and perpendicular to the center line after 
loading, and (e) that the lateral deflections are small compared to the 
length of the column. 

The present analysis Is very similar in its basic assumptions to 
that of Higgins (ref. 5 ) which also deals with the solid rectangular- 
section column; however, the two developments are entirely different, 
with the present development making evident the significant parameters. 


SYMBOLS 


A, B, E, K 
b 

Il> 12 


n 


material constants 
thickness of column 

certain integrals, across column thickness, of functions 
of (i at mldhelght 

length of column 

number of equal spaces into which cross-sectional thick- 
ness Is assumed divided 

load parameter, - 


1 - J2_ 


°E 




functions of the values of \ij_, defined by equation (22) 
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^cr 

t' 

^ cr 

w 

z 

z lA Z 2 J * 

5 0 

»1 

€ 

~ 

M- 

M-2^ * 

£ 

cr 

a 

°E 


straightness parameter. 


6Sc £0 g E 
K b crg-a 


alternate straightness parameter. 


e 


6Bct 6 0 
K h 


time 


column lifetime, value of t at which lateral deflec- 
tion become infinite 

time parameter, ^ e® 0 ’/'^ 

( AEB \ X ^ K 

lifetime parameter, t cr ) e 



width of column 

distance measured from center line of column, positive 
toward convex side 

. . z n coordinates z of centers of the n spaces (starting 
from concave side) 

amplitude of initial sinusoidal deflection 
midheight lateral deflection (not including & q ) 
compressive strain 

average compressive strain at midheight of column 

. . values of p at centers of spaces 
dummy coordinate representing z 
compressive stress 

average compressive stress on column 


plastic- strain parameter, 


f( E 


Euler buckling stress of column, — 

12(L/b) 2 
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ANALYSIS 

Assumed Law of Material Behavior 


Because the longitudinal stress in any fiber of a continuously 
deflecting column varies with time, the description of the material 
given by the assumed constant-stress creep law mentioned in the 
introduction 


e = — + Ae Bo t K (l) 

E 

is insufficient for purposes of analyzing the creep behavior of a 
column. Shanley's engineering hypotheses of creep (ref. k) can be used 
to derive the following generalization of equation (l), which covers 
behavior under varying compressive stress: (See pp. 460-A61 of ref. 1 

for the derivation.) 


e = - + f (a,e) 


(2) 


where 


f(a,s) 


K(Ae Bg ) 1/,,K 



and the dots denote differentiation with respect to time. 


( 3 ) 


Equation (2) implies that the material behaves elastically under 
an infinitely rapid increase or reduction in stress. If equation ( 2 ) is 
rewritten with the term 6 / E transposed to the left-hand side, f(cr,e) 
is seen to represent the creep rate, or rate of growth, of the plastic 

strain e - — . 

E 


Equation ( 3 ) implies that, if after a period of creep under com- 
pressive stress the stress is reduced to tension (0 negative), the 
creep rate will still remain positive - that is, compressive - although 
it may become exceedingly small as the stress takes on larger tensile 
values. The neglect of tensile creep will introduce some unconservatism 
into the analysis when the lateral deflection gets so large that the 
stress becomes tensile on the convex side of the column. The degree of 
unconservatism, however, should be small in many practical situations 
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in which the initial curvature is reasonably small and the average 
applied stress is reasonably high but not close to the Euler stress. 

In such cases the t ime between the occurrence of appreciable tension 
and the calculated ultimate collapse of the column will be a small part 
of the total column lifetime. 


Qualitative Description of Column Behavior 

The column, before application of load, is shown in figure l(a) . 

The cross section is rectangular, of thickness b and width w. A 
slight initial curvature is present in the form of ha lf a sine wave of 
amplitude 60 . 

A load, less than the Euler load, is instantaneously applied and 
produces an average stress a and an instantaneous static deflection 
6-j_(0), as shown in figure l(b). By hypothesis, the material behaves 
elastically during instantaneous stressing; therefore, the conditions 
existing immediately after load application (t = 0) are obtainable from 
an elastic analysis. At that instant no plastic strain exists anywhere 
in the column. 

With the load held constant, the column will continue to deform 
because of creep and will deflect laterally because of larger creep 
strains on the more highly stressed, or concave, side than on the convex 
side (fig. 1(c)). The problem Is to determine the history of the 
deformations . 

That the present analysis should predict a finite lifetime for 
nonvanishing initial curvature can be established by the following 
intuitive argument: A solid rectangular- section column is certainly 

less stiff than the H-section column formed from it by concentrating 
half of the material of the column at each extreme fiber, al 1 other 
conditions remaining unchanged. According to the analysis of reference 1, 
however, which is based on the same assumptions as the present analysis, 
the lateral deflection of any slightly curved H-section column approaches 
infinity in a finite time; hence, the deflection of the less stiff solid- 
section column must certainly become infinite in finite time. 


Development of Basic Equations 

For a rigorous analysis, equilibrium, compatibility, and the law 
of material behavior must be satisfied at all sections of the column at 
every instant. For simplicity, however, the deflection shape is assumed 
to be merely a magnification of the initial shape - that Is, half a sine 
wave. Because of this reduction to a single degree of freedom in the 
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deflection shape, the requirements of equilibrium, and so forth, can at 
best be satisfied only in some average way over the length of the column. 
An even simpler solution, that followed herein, is to satisfy the require- 
ments at only one section of the column, for example, at the midheight. 

Figure 2 shows the upper half of the column of figure l(c) as a free 
body. The compressive stress 0 at the cut section varies from fiber 
to fiber and also varies in time - that is, cr = a(z,t). Similarly, the 
compressive strain € = e(z,t). 


Equilibrium of vertical forces requires that . 

pb/2 


J-b/2 


cr dz = ab 


00 


Equilibrium of moments requires that 

~b/2 
-b/2 


Pb/2 

crz dz = -cb (S 0 + Bt) 

J-b/2. 


(5) 


The assumed law of material behavior (eqs. (2) and (3)) gives the fol- 
lowing equation for the strain rate : 


e = 


a K\Ae 
E 


U Bg ) 


l/K 


0 

6 " E 


I ' 1 


( 6 ) 


Finally, the assumptions that the shape is sinusoidal, that plane sec- 
tions remain plane and perpendicular to the center line of the column, 
and that the deflections are small result in the following relationship 
between the strains and the lateral deflection at the column midheight: 




e = -z 


1 + e 


120tTi6*| _ 

-Z + € 

Eb 2 


(T) 
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wiiere e is the average strain across the section, or 


e = ( e )z=0 

and cfg is the Euler huckling stress of the col umn ; that is. 


12(L/b) 2 


(8) 


(9) 


Reduction of Basic Equations 

Equations (4) to (j) can he reduced to a single equation in terms 

of the elastic strain e - — at the midheight of the column, which is 

E 

also a function of position and time . Equations (if), ( 5 ), and (6) are 
first rewritten as follows (with cognizance being taken of eq. (7)): 
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b_( a) = KA 1 / K e Bo ’/^- 

BtV " eJ i_; 



Equation (12) can be reduced further if e is eliminated through. use 
of equation (10) and S-j_ eliminated through use of equation (ll). The 
result, after slight simplification, is the following equation in which 

e - ^ is the only unknown: 
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where i is a dummy variable representing z. Introducing the following 
parameters : 

the plastic -strain parameter 



the time parameter 



the initial- straightness parameter 

6Ba B 0 °E 

“ K h ap-a o 

e ^ = S 


and the load parameter 



permits equation (13) to he 


simplified to 


0 *) 


(15) 


(16) 


(17) 


J- 1/K = s 2 Z /b e 
St' ^ 



p 4 + 


12zP 


Pl/2 




(18) 


where p = p(z,t) is subject to the initial condition p(z, 0) = 0. 

Equation (l8) gives, essentially, the rate of growth of the plastic 
strain to the l/K power for any fiber of the column as a function of 
the position z of that fiber, the instantaneous value of the plastic- 
strain parameter p in that fiber, and the two integrals involving the 
instantaneous values of p in all the fibers of the column. 
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For the initially perfectly straight column (S = 1 ), equation (l8) 

will he satisfied by p = t'^, as is expected, inasmuch as p = t'^ is 
equivalent to equation (l), the law of material behavior for direct 
compression. 


Numerical Method of Solution of Equation (18) 

Equation (18) lends itself readily to numerical solution by means 
of the modified Euler method (ref. 6). The following adaptation of this 
method has been found suitable for use in conjunction with a high-speed 
digital computer, the National Bureau of Standards Eastern Automatic 
Computer (SEAC): 


Let the cross section of the column at the midheight be divided into 
n equal spaces. The width of each space will be b/n and the centers 
of the spaces, starting at the left-hand, or concave, side of the column 


and going toward the right, will be located at z =* 



1 b 

2 n' 


z , z = t + 5 b 
2 2 2 n' 


Let the values of the plastic-strain param- 


eter p at the centers of these spaces be denoted by p-^, p2, • • • P n , 


respectively. The two integrals appearing in equation (l8) can then be 
approximated in terms of the parameters p-j_ and z^, the very simplest 


approximation for machine -calculat ion purposes being that obtained by 
assuming that the integrand is constant over the width of a space and 
equal to its value at the center of the space. Then, 




J(p x + p 2 + 



( 19 ) 


and 



( 20 ) 
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Writing equation (18) for the specific z = z ± , z 2 , • • • z n gives 

the following system of ordinary differential equations which is approxi- 
mately equi valent to equation (l8) : 


d 



dt' 







dt' 



n e L 


-^2 + Il+6(-l+|)pi2 


r (2D 



-1+ 


2n-l 
n _ 


- fe +I 1+ 6(-l+ ^i)PXa 


Letting 


= S 


-1 + 


21-1 

B - 


-Hi+I^^-l+^^PIg 


simplifies equations (21) to: 


d (^ 1 1 / K 

dt' 


= R, 


d [ q 2 l/K ) 

IP ^ 


4n l/K ) 


dt 1 


= R 


i n 


( 22 ) 


f (23) 
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The parameters Rq are functions of t' "because they involve the 

Hq(t'). 

Equations (23) define curves of the pq plotted against t' and 
can "be solved in a step-by-step fashion subject to the initial conditions 
p i (0) =0. As may be expected on the basis of the physical argument 

presented previously, at least one of the pq, probably pq, will tend to 
approach infinity asymptotically at some finite value of t'. A step- 
by-step solution of equations (23) on the basis of equal increments 
of t' therefore should not be attempted. Rather, the solutions should 

be carried forward on the basis of equal increments in pq or pq^/^, 

and the increment in t' should be regarded as one of the unknowns in 
each step. A detailed account of the calculation procedure follows. 

Suppose that the curves of pq, p£, * • • Ph plotted against t' 
have been traced up to some particular value of t' . It is now desired 
to carry these curves forward one more step corresponding to a specified 
finite increment in pq denoted by Apq. To be determined are the 
increment At' in t' and the increments Ap^ Apj , • • • Ap n ’k* le 

other pq corresponding to the specified increment Apq in pq. Since 

the values p-|_ and Apq are known, the increment A^pq 1 / K j in the 

quantity p-^/K can be computed. The first of equations (23) then 
gives the following approximation to At': 


At' 



(2k) 


With At ’ 

a(, 3 1/K > • 


known, the remaining equations give the increments 
• • aLWk as follows: 



A(p 2 1 / K ) = At’RgCt') 

A (p 5 1/,K ) = Afc'R 5 (t') 

• • * • • 

A (^n l/K ) = At,R n(t') 




(25) 
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From these increments, the values of p 2 , p^, • • ■ Mn corresponding 

to t' + At' can he computed. The curves have now been tentatively 
advanced one step and tentative values of Rq(t' + At') can he deter- 
mined. Improved values of At' and A^p 2 ^/^j , A^p^/^j, • • • A^p n ^^ 

can now he obtained by replacing the right-hand sides of equations (24) 
and (25) by averages based on the values of Rj_ at the end of the 
interval as well as at the beginning of the interval, that is, by using 
the equations 


At' 


(u 

1 1 

-j- 

\ 1 ) 2 

R 1 (t' ) ' R x (t ' + At' ) 


(26) 





At' 


1 

2 


R 2 (t’ ) + R^t' + At' ) 



At' 


1 

2 


R 5 (t' ) + R 5 (t' + At' )! 


r (27) 



At' 


1 

2 


R n (t ' ) + R n (t ' 


+ At' ) 




The improved values of Pp(t' + At') obtained from equations (27) 

give new values of R^(t' + At') which can be substituted in equations (26) 
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and (27) in order to obtain still better values of At' and 

This process can be repeated until no change 
is evident in the values of At' and Pg(t' + At'), p^(t ' + At'), • • • 
^(t* + At'). 


The procedure just described can also be used without modification 
to start the curves. It is necessary to use only the fact that m(0) = 0. 


Approx ima te Solution of Equation (l8) Based on n = 2 


If, for purposes of numerical integration, the cross section were 
divided into two equal spaces only, then the system of equations (21) 
would reduce to the following pair of equations: 



dt' 


d (^ l/K ) . l) 

dt ' 




> (28) 


■J 


With the differentiations indicated on the left-hand sides carried out, 
these equations become identical in form with the equations which arise 
in the analysis of the two-flange column - namely, equations (21) and (22) 
of reference 1. The solution presented in reference 1 can therefore be 
used to obtain the approximate solution of the present problem corre- 
sponding to n = 2. In order to use the solution in reference 1, certain 
substitutions should be made in the symbols and graph labels of that 
paper and, then, the results of that paper will apply directly to the 
present problem. The substitutions are as follows: 

(1) Replace 7 R /7 L toy S 

(2) Replace q L by Pqp| - 
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(3) Replace t] r "by P 2 (^ " |) 

(4) Replace 7 L t by ^S -1 / 2 ! 1 

(5) Replace 7 L t cr by (j^ - i) 1 ^' 1 / 2 ! ' cr 

(6) The constant K is unchanged 

Lateral Deflections, Strains, and Stresses 

With the pj_(t') completely deter min ed hy the solution of equa- 
tion (l8), the time history of the plastic strains e - ^ is, in effect, 

known. Equation (ll) can then he used to compute the midheight lateral 
deflections S-,(t), and equation (10) gives the history of the average 
midheight strain if(t). The strains e(z,t) are then obtainable from 

equation ( 7 ). A knowledge of e - ^ and e is sufficient to determine 

the midheight stresses cr(z,t). 

RESULTS AND DISCUSSION 
Significant Parameters 


The usefu ln ess of data obtained from tests is enhanced if the data 
are presented in terms of the appropriate parameters. Equation (l8) of 
the present analysis suggests that test data on the creep life of columns 
whose material creep curves are adequately represented by equation (l) 
may profitably be plotted in terms of the parameters K, t' cr , S, 

and o/cr-g, where t' cr is the value of t 1 corresponding to collapse of 

the column. If the analysis were perfectly applicable, the test data 
would give a single curve of t' cr plotted against S for a given value 
of ayop; and a given value of K. Scatter of test points away from such 
a curve will, of course, occur, because of shortcomings of the theory - 
for example, the neglect of third-stage creep in equation (l), the approxi- 
mate engineering hypotheses used to generalize equation (l), the neglect 
of tensile creep in the generalization, the assumption of an exactly 
sinusoidal shape for both the initial and the subsequent deflections, 
and the assumption of instantaneous elastic loading. 
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Inasmuch as the parameter S contains the parameter cr/ag and each 
plotted curve or test point has a definite value of ct/oe associated with 
it, a simpler parameter S' may be used instead of S, where 

6 B 0 £0 
S' = e K 


1 - 


= S 


°E 


(29) 


If the test results for a single material at a single temperature 
are to be plotted, then A, B, E, and K are the same for all the 
tests and the significant parameters may be simplified further to 


t cr e Ba / K , a and 0 / 0 %. 


Calculated Values of t' cr 

In the section entitled "Analysis" a numerical method of solving 
equation (l 8 ) was given, based on dividing the column cross section at 
the midheight into n equal spaces and writing equation ( 18 ) for the 
center of each space. This method has been carried out and found to be 
satisfactory for the case K = 1, oya-g = 0.7 and 0.9; and for several 

values of S. The value n = 14 was found to lead to sufficiently 
accurate results in a trial calculation, and this value was used for 
all the cases. 

Figure 3 shows a typical curve of pp plotted against t'. (The 
curves of pg* M-j; * • • are not shown, but they would lie suc- 

cessively below the curve of pp.) This figure tends to confirm the 
intuitive argument presented previously to the effect that the deflec- 
tions should approach Infinity in a finite time. Computation was stopped 
when the computed increment in t' corresponding to the specified 
increment in pp became less than the precision of the computation. 

The last computed value of t' was taken as being equal to the asymptotic 
value t' cr for all practical purposes. 

The values of t' cr taken from a number of curves like the one in 
figure 3 are summarized in figure 4 in plots of t' cr against S for 
ofa Q = 0.7 and 0.9. The circles are the actually computed data through 
which the curves are faired. 

As is to be expected, t' cr tends toward zero for columns that are 
more and more initially curved (S — > 0) and toward infinity for columns 
approaching straightness (S — >•!). In order to show more clearly the 
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"behavior of the curves in these two regions, the data of figure k- are 
replotted in figures 5 and 6 in which logarithmic scales are used to 
magnify the region of S = 0 and S = 1, respectively. 

All the computations were performed on the National Bureau of 
Standards Eastern Automatic Computer. The computation of a value of 
t' cr took about 15 minutes, exclusive of the exploratory computation 
time required to establish optimum values of precision, the step size 

and- the number of cross-sectional spaces n. The calculated values 

of t'cr are felt to be accurate to within 2 percent. 

An approximate solution for the solid rectangular- section column 
corresponding to n = 2, as is indicated in the section entitled 
"Analysis," already exists i mp licitly in the two -flange -column analysis 
in reference 1. Figure 5 of reference 1, which summarizes the theoret- 
ical lifetime data for two-flange columns, can therefore be used to yield 
approximate information on the lifetimes of solid columns by replacing 
the label 7^/7p "by the label S and the label 7 pt cr by the label 



l/K -1 / 2 
S t 


cr 


In figure J, the approximate results thus obtained for the 
rectangular- section column for the cases where K = 1 and a/cr-g =0.7 
and 0.9 are compared with the much more accurate results for these cases 
based on computations with n = ik and obtained from figure 4 of the 
present paper. 


CONCLUDING REMARKS 


On the basis of certain assumptions regarding the shape of a column 
and the creep behavior of the material under varying stress, an analysis 
has been made of a slightly curved solid rectangular- section pin-ended 
column carrying a constant load. The analysis led to a differential 
equation for the plastic strains at the midheight cross section as a 
function of time. The form of the equation indicates the significant 
parameters of the problem, which parameters may be useful in plotting 
test data on the creep life of columns. These are a lifetime parameter 
t' cr , an initial- straightness parameter S or S', and the ratio of the 
average applied stress to the Euler stress cr/cr-g. 

A n ume rical method of solution of the differential equation, suitable 
for use with a high-speed digital computer, has been described, and • 
typical computed results have been presented. 
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The existence of a finite lifetime, although not evident from the 
differential equation, was argued intuitively and confirmed by the 
numerical computations. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 31, 1953- 
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(a) Before loading. 


(b) Immediately after 
application of 
load (t-O). 


(c) At time t after 
application of 
load. 


Figure 1.- Schematic description of column. 



abw 



6Q + fil(t) 


Figure 2.- Upper half of column of figure l(c) . 
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Figure 5 


t' 


, - Variation of plastic-strain par- 


parameter t'. K = lj 5 /ct e = 0.7; 










Figure 7 .- Comparison of approximate solution, based on n = 2, with 
accurate solution, based on n = Ik, for K = 1 . 
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